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ABSTRACT 

Three  K(9)  functions  were  used  to  describe  hydraulic  conductivity 
data  from  a  layered  field  soil.  The  functions  were  KQ(0/Qm)&, 
K0((e-ec)/(em-ec))n  and  K0exp{a(e-em)}  where  K0  =  K0(z)  =  g(z)Km, 
g(z)  is  a  scaling  factor  that  varies  with  depth  z,  Km  a  constant 
and  9CI  9m,  P,  n,  and  a  are  parameters.  For  each  function  three 
cases  were  considered:  Case  1  fit  discrete  values  by  depth  to 
g(z)  and  a,  p  and  n:  The  second  case  treated  the  scaling  and 
exponential  parameters  as  continous  functions  of  depth:  And  the 
third  case  fit  scaling  factor  as  a  discrete  function  of  depth  but 
held  p,  n  and  a  constant.  More  variation  in  r2  and  MSE  was  found 
between      cases      than      between      functions. 


INTRODUCTION 

Most  soil  physics  field  studies  to  date  have  been  directed  at 
obtaining  detailed  descriptions  of  the  various  soil  horizons  or 
alternatively  the  course  of  action  has  been  to  ignore  layering 
and  treat  the  soil  as  being  a  uniform  body.  Describing  a  soil 
horizon  by  horizon  is  a  tedious  task  with  no  clearly  defined 
stopping  point.  Warrick  et  al.  (1977)  used  scaling  to  reduce  the 
amount  of  effort  required  to  describe  the  hydraulic  properties  of 
small  laboratory  cores.  The  objective  of  this  paper  is  to 
evaluate  scaling  as  a  technique  for  describing  hydraulic 
properties      of      a      layered      field      soil. 

THEORY 

Spatially-varying  hydraulic       conductivities       can       be       approxinated       as 

the      product      of  a      function      of      depth      and      a      function      of      volumetric 

water       content  (Warrick       et       al.,       1977),       i.e.: 

K(9,z)     =    gi<  (9)  [1] 

where  K(9,z)  is  the  spatially-varying  hydraulic  conductivity, 
g=g(z)  is  the  scaling  factor  expressed  as  a  function  of  depth 
(z),  and  K(9)  is  the  nominal  hydraulic  conductivity  and 
depends  only  on  water  content.  Taking  the  log  of  both  sides  of 
Eq.       [1]       results       in 

iog{K(9,z)]    =    log[g(z)J    +    log{K(9)]  [2] 


Choose 


G       B 
K(9)     =     Km(  — )  [3] 


e  -ec  n 

K(9)     =    Km( )  [4] 

K(9)     =    Km    exp{a(9-em) 1  [5] 

where      Km      is      the      hydraulic      conductivity      at      the      reference      water 
content      0m.      ec,      B,      n,      and      a      are      parameters.      Eq.      [3],      [4]      and      [5] 
were       taKen       from       Watson       (1967),       Brooks       and       Corey       (1964)       and 
Davidson        et        al.        (1963),        respectively.        Substitution        of        one        of 
Eq.      [3],      [4],      or      [5]      into      Eq.      [2]      results      in 

9 
log[K(z,e)]    =    log{K0(z)J    +    3-log( — )  [6] 

em 

©    — © 
log{K(z,©))    =    log{K0(z)J    +    n • 1 og ( ]        [7] 

em"ec 
ln[K(z,e)]     =     ln[K0(z)]     +    a(6    -6m)  [8] 

where      K0(z)      =      g(z)Km. 

While  no  formal      statistical      test      is      known      for      choosing      the      best 

g(z)  or  functional       form       of       K,       a       qualitative       appraisal       can       be 

made  on  the     basis     of     r2     and     mean     square     errors.     Three     cases     for 

each  of  Eq.      [6],      [7],      and      [8]      will      be      considered. 


CASE  1: 

log{K0(z) ]  =  log{K0(z, ) ) 

Pj=  P(Zj);  n,=  n(2j);  and  a,=  a(z,) 

where  Z\     denotes  a  depth  where  measurements  were  made  and  g, ,  0j  , 

a,    and   n,  are   discrete  estimates.   This  case   required   fitting 

functions  at  each  individual  depths. 

CASE  2: 

log{K0(z)  )  =  log[K0(z)  ) 

P  =  P(z);  n  =  n(z);  and  a  =  a(z) 

where    log[K0(z)j    =    a0    +    a(z    +    a2z2    +a3z3    +a4z4    +    a5z5 

cc(z),3(z)    or    n(z)     =    b0    +    b|Z    +    b2z2    +  b3z3    +b^z4    +    b5z5 

and  a0 , a 1 a5  and  bg,b| , . . . ,bg  are  regression  coefficients. 

CASE  3: 

log(K0(z)  i  =  log[K0(Z|  )  J 

P,  n  and   a  are  constant  over  al I  depths. 

Case  3   results  in   a  single   relation  between   K  and   0,  K (©) , 
which  will   be  hereafter   referred  to   as  the   nominal   hydraul  ic 
conductivity  function. 

Mater  i  a  1 s  and  Methods 

The  field  study  was  carried  out  on  Mu i r  silt  loam  (fine-  silty, 
mixed,  mes i c  Pachic  Haplustoll)  located  on  Kansas  State 
University  Ashland  Research  Farm  about  10  km  south  of  Manhattan, 
Kansas.  The  Mu i r  series  consists  of  deep  nearly  level  soils  on 
river  and  creek  terraces  formed  in  deep  alluvium.  The  3  by  4  m 
test  plot  was  bermed  with  soil  and  sand  bags.  An  aluminum 
access  tube,  4.13  cm  outside  diameter  and  183  cm  long,  was 
installed  to   a  depth   of  168   cm  in   the  center   of  the    plot  4 


months  before  starting  infiltration  on  14  July  1983.  Mercury- 
manometer  tensiometers  were  arranged  in  a  circular  pattern,  80 
cm  in  diameter  around  the  access  tube  at  the  0,  20,  40,  60  ,  80, 
100,  120,  140,  and  160  cm  depths.  Water  was  maintained  at  a 
depth  of  5  to  6  cm  on  the  plot  until  soil  water  pressure  and  © 
appeared  constant  with  time.  The  plot  was  then  covered  with 
plastic  and  2  to  3  cm  of  soil  to  prevent    evaporation. 

Neutron  moisture  observations  were  made  at  12  depths  10,  20, 
30,  40,  50,  60,  70,  80,  100,  120,  140,  and  160  cm  as  the  mean  of 
two  64s  counts.  During  the  early  drainage  phase,  readings  were 
taken  every  2  to  4h  for  the  first  12h  then  less  frequently  at 
larger  times.  Soil  water  content  measurements  were  made 
immediately  after  the  tensiometer  readings.  Hydraulic 
conductivities  were  estimated  using  a  modified  instantaneous 
profile  method  described  by  Rose  et  al.  (1965). 

Percentages  of  sand,  silt,  and  clay  were  determined  using  the 
hydrometer  method  (Day,  1965).  Except  for  the  upper  15  cm  and 
the  lower  30  cm  (135  -165  cm)  where  the  soil  texture  is  silt 
loam  the  soil  profile  is  predominantly  a  silty  clay  loam.  Average 
particle-size  distribution  for  each  depth,  given  in  Table  1, 
shows  that  the  mean  fraction  of  clay,  silt,  and  sand  were  28, 
54,  and  18/,  respectively.  Table  1  also  shows  that  the  average 
bulk  density  value  for  the  top  100  cm  is  1.4  Mg-m"3;  this  value 
decreased  to  a  minimum  at  140  cm  with  the  lower  60  cm  of  the 
profile  having  an  average  density  of  1.32  Mg-m"3. 


The  neutron  probe  was  calibrated  by  regressing  volumetric  water 
content  vs.  neutron  probe  count  using  data  from  borings 
obtained  from  the  plot  area  following  the  infiltration-drainage 
experiment.  To  extend  the  range  of  the  field  calibration  to 
the  wet  end  the  area  was  reflooded  and  soil  samples  taken.  Bulk 
densities  and  volumetric  water  content  were  estimated  from  7.6  by 
7.6  cm  undisturbed  cores  taken  with  a  thin-walled  hydraulic 
probe.  There  were  a  few  exceptions  where  the  average  bulk  density 
was  used  to  convert  gravimetric  water  to  a  volume  basis.  The 
narrow  range  of  the  clay  content  and  bulk  density  (Table  1) 
suggested  using  one  neutron  probe  calibration  curve  for  all 
depths  (Fig.  1).  Observations  at  10  cm  depth,  where  the 
calibration  was  affected  by  the  restricted  soil  volume  near  the 
surface,  were  deleted  from  analysis. 

Chemical  analysis   of  the   well  water  used  in  this  study  was    as 
follows:   0.6  mmol   L*'  Na+,   0.1  mmo  1   L"'  K  +  ,   15.2  mmo  I  L"f  Ca2  +  , 
3.6   mmol  L"1    Mg2+,  0.21  mmol  L'1   CI",  1.4  mmol  L_l   S042-,     0.2 
mmol    L"'   N03",  10  mmol  L"1  HCO3",  pH  =  7.08,  and  EC  =  0 .  85  os -m"1* 

Resu 1 ts  and  D  i  scuss  i  on 

Water  content  profiles  at  0.0,  0.2,  0.5,  1.4,  3.4,  7,12,  18. 3, and 
28. 4d  are   presented  in   Fig.   2.   Measured   water   contents   are 
indicated  by      the  solid   circles.  Average  water  content  change 
above     70   cm  was   0.03  m3-m"3  during  drainage.  A  larger  change, 
0.06  m3-m"3,   occured  below   70  cm.   The    rate  of  change  of  water 
content,  ae/at,   during  early  drainage    phase  was  about  80  times 


greater  than  the  rate  of  change  of  water  during  late  drainage 
phase . 

Figure  3  shows  total  hydraulic  head  plotted  against  depth. 
During  the  28-day  drainage  period,  the  pressure  head  at 
deeper  depths,  100,  120,  140,  and  160  cm,  changed  as  much  as 
twice  that  of  the  upper  depths.  Hydraulic  head  gradients 
approached  a  magnitude  of  3  near  the  end  of  the  drainage  phase. 

Three  K (9)  were  fitted  at  each  depth  for  Case  1.  9C  i n  the 
Brooks  and  Corey  Model  was  estimated  from  curve  fitting  total 
water  above  z  (W)  data, 


W  =  W(z, t)  = 


9dz 


[9] 


0 


where      t      is      time      in      days.        For      the      Brooks      and      Corey      Model      and 

assuming      a      unit      gradient      Sisson      et      al.      (1980)      expressed      W      as, 

z  1/(n-1) 

W    =    W(z,t)     =    9cz+(1-1/n)z(9m-9c)  ( )  [10] 

At 

here  9m  was  estimated  as  the  average  of  9(z,t)  after  16  days  of 
ponding.  Note  that  the  n  here  is  1/n  in  Sisson  et  al.,  (1980).  A 
is  given  by  A  =  Km/[(9m-9c)/n],  Km  was  estimated  from  the  final 
infiltration  rate  at  the  soil  surface  after  16  days  of  ponding, 
and  9C  and  n  were  estimated  by  non-linear  least  squares  fitting 
to  Eq.  [10]  (PROC  NLIN  procedure  provided  by  Helwig  and  Council, 
1979). 
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The  three  K(0)  functions  (Eq.  [6],  [7]  and  [8])  were  fitted  for 
Case  1.  The  results  for  the  140  cm  depth  are  shown  in  Fig.  4  and 
were       considered       typical       of       most       depths.  Estimates       of       the 

parameters  are  given  in  Table  2.  The  three  functions  performed 
equally  well  in  that  similar  r2  values  were  obtained  at  each 
depth.  The  results  from  other  depths  tended  to  follow  the 
sigmoidal  shape  shown  in  Fig.  4.  None  of  the  K(9)  functions 
cosidered  here  could  mimic  such  a  shape.  The  eleven  regression 
lines  resulting  from  fitting  Eq.  [8]  for  Case  1  are  shown  in 
Fig.   5. 

The  slope  of  each  regression  curve  obtained  from  fitting  the 
three  different  K(9)  functions  for  individual  depths,  was  tested 
against  all  other  slopes  and  the  results  of  this  test  are 
presented  in  Table  3.  This  test  was  a  t-test  dependent  on  the 
difference  between  two  slopes  as  well  as  individual  standard 
errors.  This  test  gave  the  same  ranking  of  exponential  parameters 
(slopes)  regardless  of  the  K(0)  function.  Slopes  from  the  120 
and  60  cm  depths  had  the  largest  differences  (Table  2).  The 
ranking  of  exponents  remain  the  same  regardless  of  the  K(9) 
function.  It  may  be  concluded  from  Table  3  that  significant 
differences  existed  among  exponential  parameters  at  the  one 
percent        probability        level. 

Case      2      fit      the      slopes      of      Eq.      [6],      [7]      and      [8]      as  polynomials      of 

z.          The       maximum       r2       improvement       technique       was  used       to       fit 

polynomial       models       to       the       field       data       (Helwig       and  Council,       1979). 

Estimates      of      the      polynomial      coeffients      are      given      in  Table      4.      The 


continuous  polynomial       functions       produced       higher       MSE       and       lower 

rJ       values  than       Case       1       (Table       5).          The       continuous       polynomials 

offer       the  advantage      of      estimating      hydraulic      conductivity      at      any 

depth     over  the     0     -  1 60    cm    range. 

Case  3  requires  fitting  curves  of  the  same  slope  to  Eq.  [6],  [7] 
and  [8].  The  fitting  was  done  using  the  method  of  dummy 
variables  (Draper  and  Smith,  1981).  The  results  of  this  fitting 
are  shown  in  Fig.  6  for  the  Davidson  model.  A  numerical 
comparison  of  cases  is  presented  in  Table  5.  MSE  values  of  the 
single  slope  model  were  about  twice  that  of  the  variable  slope 
models.  Cases  1  and  2  also  produced  rJs  that  were  9  to  15X 
higher  than  the  single  slope  model.  The  parallel  curve  models 
explained  more  than  74/  of  the  variation  of  log{K}  or  ln{K) 
around  the  means  (Table  5).  The  basic  advantage  of  the  parallel 
model  is  that  only  one  function  needs  to  be  estimated  to 
describe       a       soil. 

Case  3  allows  a  single  K(9)  for  the  whole  profile,  once  the 
scaling  factor  g(z)  has  been  estimated.  Assuming  the  Davidson 
Model  Km  for  the  nominal  hydraulic  conductivity  was  estimated  as 
the  log  mean  of  K0(z)  over  all  depths  and  is  shown  in  Fig.  7  as 
the  solid  curve.  Hydraulic  conductivity  data  from  all  depths  were 
adjusted  by  the  difference  between  lnK0(z)  and  lnKm  and  also 
plotted  on  Fig.  7.  A  regression  of  the  data  in  Fig.  7  indicated 
an  r2  of  0.91+  considereably  higher  than  most  of  the  r2s  given  in 
Table     2     and     any     of     the     r2s     in     Table     5.     If     the     r2s     in     Table     5     had 
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not  been  adjusted  for  the  mean  of  the  dependent  variable  they 
would      also      exceed      0.90. 

Further  field  work  would  be  required  to  establish  how  reliably 
the  single  curve  predicts  hydraulic  properties  of  the  Muir  silt 
loam         other         locations. 

Summary         and         Conclusions 

1.        Three        K(9)        were        tested        their        ability        to        describe  field 

measured      K(9)      data.      Exponents  in      the      three      K(9)      functions  were 

allowed      to      vary      with      depth      in      two      tests      and      were      fixed  to      a 

single      value      in      a      third      test.  All      three      K(9)      functions      performed 

equally        well        as        determined  by        r2,        although        the        r2  for 

the      variable      and      fixed      models  were      different. 


2.  When  the  exponents  in  the  three  K(6)  functions  were  fitted 
as  discrete  and  polynomial  functions  of  depth,  the  discrete 
functions  yielded  higher  r2  (0.89).  When  one  exponent  was 
used  for  all  depths,  the  r2  ranged  from  0.74  -  0.78.  The 
advantage  of  the  polynomial  procedure  was  at  points  between  the 
depths  studied  could  be  interpolated.  The  advantage  of  fixed 
slope  model  is  that  only  one  exponent  needed  to  be  determined 
for  all  eleven  depths.  The  discrete  model  had  the  highest 
precision. 
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Fig.  2.  Soil  water  content  profiles  of  Hulr  silt  loam 
during  drainage  cycle. 
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Fig.  4.  Hydraulic  conductivity  volumetric  Mater  content 
relatione  of  the  140  cm  depth.  Solid  lines  were  obtained 
from  fitting  the  discrete  slope  and  Intercept  functiona 
to  field  measured  K-0  data  and  the  circles  refer  to 
instantaneous  profile  method. 
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Fig.  5.  Hydraulic  conductivity  volumetric  water  content 
relatione  for  all  depths.  Solid  lines  were  obteined 
from  fitting  Eq.  [8]  ee  e  discrete  model  to  the  field 
measured  K-0  deta. 
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Fig.  6.  Hydraulic  conductivity  volumetric  water  content 
relatione  for  all  depths.  Solid  lines  were  obtained 
from  fitting  Eq.  [B]  as  a  fixed  exponent  model  to  field 
measured  K-0  data. 


25 


rT 


o 
o 

> 


EEEEEEEEUUU 
uuuuuuuu 

ooo 

OOOOOOOOOfllt 


-H(U(n^inio^aioio< 


iim  ■  i  ■ 

M 
O 


■in  iii    i 


|iui  i 

o 
O 


inn  1 1  i 


m 

E 


,18  * 

z 
o 
u 

cr 

LU 

I- 
< 


-m 


u 

M 
UJ 


C\J 


ni 

i 
O 


in 

i 
O 


(P/UI3)    AiiMiDnaNOD  DiinvaaAH 


Fig.  7.  Noalnal  hydraulic  conductivity  based  on  Eq.  [8]  . 
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mean  aquare  error  (MSE)  obtained  from  fitting  the  three 
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Literature        Review 

Understanding  of  the  water  movement  into  and  through  soil 
profile  is  of  significant  importance  because  it  is  the  basic 
input  for  irrigation  and  drainage  systems  design,  soil  management 
practices,  crop  selection  and  land  use,  and  environmental 
aspects.  These  practices  are  all,  to  a  certain  extent,  dependent 
upon  the  hydraulic  properties  of  soil  which  include  the  water  and 
conductivity  chracteristic  that  relates  soil  water  content  with 
pressure  potential  and  hydraulic  conductivity,  respectively. 
Water  infiltration  and  redistribution  are  time  dependent  and 
essential  for  solute  transport  through  soil.  Gardner  and  Widtsoe 
(1921)  developed  an  equation  to  relate  the  mean  values  of  water 
content  following  irrigation  with  time,  T,  in  days 
W     =     14.6     +     7.6     e-<MKT     +     7.6     e"0^1 

where  W  represents  the  mean  value  of  water  contents  measured  to  a 
depth  of  183  cm  and  expressed  as  a  percentage  on  a  dry  weight 
basis.  While  considering  the  physical  processes  involved  in  loss 
of  water  from  a  field  plot  of  a  uniform  sandy  loam  after 
irrigation,  Richards  et  al.  (1956)  found  that  the  total  water,  W, 
above  a  given  depth  could  be  closely  related  to  time  by  an 
equation  of  the  form 
W  =  aT-0 

where      a      and      b      are      costants,      and      that      the      rate      of      loss      of      soil 
water       was      inversely       proportional       to       time,       T,      i.e; 
dW/dT   =    -abT^1  =      -bW/T 
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Using  the  same  equation,  Wilcox  (1959)  noted  that  straight  lines 
were  obtained  when  log  W  was  plotte  i  vs.  og  T.  Also,  he  reported 
that  the  fit  was  quite  close  except  during  the  early  drainage 
phase  and  that  in  any  one  soil  equal  water  contents  are 
accompanied  by  equal  rates  of  drainage.  Richards  and  WeeKs  (1953) 
suggested  a  method  to  calculate  the  hydraulic  conductivity  from 
data  obtained  during  transient  changes  of  water  content  and 
tension  in  soil  column.  The  basis  of  their  calculation  was  the 
Darcy       equation 

q=Ki 
where       q  =  flux,       K       is       the       hydraulic       conductivity,       and       i       is       the 
hydraulic      head      gradient.      They      rewrote      the      equation      in      the      form 

1/A      *      dQ/dT      =      Ki 
where      dQ/dT      is      the      rate      of      flow      of      water      past      a      given      cross 
section       area,       A,       of       the       soil       column.       The       derivative       dQ/dT       was 
evaluated        using        a        method        suggested        by        Richards        (1938)        to 
calculate      the      volume      of      water,      Q,      in      a      soil      column      as      follows: 


Q    =    A 


z 

P  gdz 

O    b 


where  P^  is  the  soil  bulk  density,  A  and  z  are  the  cross 
sectional  area  and  length  of  the  soil  column,  respectively,  and 
g  is  the  gravimetric  water  content.  This  integral  was  evaluated 
at  different  times  during  desorption  along  the  soil  column  to 
calculate  the  flux  (dQ/dT  =  q),  then  the  flux  was  divided  by  the 
hydraulic  gradient,  obtained  from  tensiometers,  to  calculate  the 
hydraulic  conductivity.  The  method  of  Richards  and  Weeks  was 
modified      by      Richards      et      al.      (1956)      and      Ogata      and      Richards      (1957) 
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for  the  analysis  of  field  observations  to  determine  the  relation 
of  water  content  to  suction  and  hydraulic  conductivity  of  a  fine 
sandy        loam        of        uniform        profile.  However,        they        expressed 

reservations  on  the  possibilty  of  in  situ  measurement  of  the 
hydraulic  conductivity  of  unsaturated  soil  profile  that  was  not 
uniform. 

During  the  late  50Js  the  most  reliable  work,  concerning  the 
rate  of  drainage  following  irrigation,  was  that  in  which  plots  of 
bare  soils  were  covered  to  prevent  evaporation  and  to  prevent 
rainfall         from         falling,         due         to  the         unavailability         of         an 

independent       measurement        of        evaporation.  Wilcox       (1959)       studied 

this  phenomena  on  bare  soils  which  were  irrigated  with  sufficient 
water  to  wet  the  soils  to  field  capacity  and  then  covered  to 
prevent        evaporation.  The        results        indicated        that        the        finer 

textured  the  soil  the  greater  were  the  rates  of  moisture  loss, 
however,  slope  of  the  line  which  relates  the  rate  of  moisture 
loss  to  time  was  much  less  with  a  clayey  than  with  a  sandy  soil. 
He  also  noted  that  the  moisture  content  following  desorption 
increased  progressively  with  depth  because  with  increasing  depth 
each  successive  layer  of  soil  received  more  and  more  water  from 
above. 

While  tensiometers  are  the  most  common  and  satisfactory 
devices  to  measure  soil  water  tension  and  hence  the  hydraulic 
head  gradient,  the  neutron  method  of  measuring  soil  water  has  the 
advantages  of  precision  and  cost  compared  to  gravimetric  sampling 
because       it       satisfies       the       most       common       requirement,       the       non- 
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destructive      attribute,       by      allowing      repeated      sampling      in      the      same 
access       hole       when       measuring       changes       in       stored       water.  This 

technique      was      used      by      Burrows      and      Kirkham      (1958)      and      Nielsen      et 
al.      (1959)      to      measure      changes      in      stored      water      of      a      soil      profile. 

Rose  et  al.  (1965)  represented  a  theory  based  on  the  water 
conservation  equation  for  a  given  volume  of  vegetation-free  soil 
to  determine  the  hydraulic  conductivity  in  the  field  as  a 
function  of  depth  over  the  entire  range  of  water  contents  on  a 
soil       of        non-uniform        profile.  This        method        is        known        in        the 

literature  as  the  instantaneous  profile  method.  On  using  this 
method,  there  are  three  options  with  respect  to  the  measurement 
of  soil  water  content  and  soil  water  pressure  head  profiles:  i) 
in  situ  measurement  of  water  content  and  pressure  distributions, 
ii)  in  situ  measurement  of  water  content  distribution  and 
inferred  pressure  head  from  water  retention  data  ,  and,  iii)  in 
situ  measurement  of  pressure  head  and  inferred  water  content 
from  retention  data.  The  first  option  is  inherently  the  best 
choice  Klute  (1972),  because  the  in  situ  determined  0(h)  curves 
often  disagree  with  8(h)  curve  determined  on  undisturbed  core 
samples  collected  from  the  same  site.  This  disagreement  has  been 
more  experienced  with  fine-textured  soils  (Luxmoore  et  al. 
1981),  but  it  also  occurs  for  the  coarse-textured  soils  (Dane, 
1980).  At  one  location,  reasonably  accurate  and  precise  values 
can  be  obtained  after  infiltration  during  the  redistribution  and 
drainage  of  water  within  and  from  soil  profile  where  measurements 
were  made  for  both  soil  water  contents  and  soil  water  pressure 
(Nielsen      et      al.        1964       ;       ran      Bavel      et      al.        1968a). 
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Wilcox  (1960)  investigated  the  effect  of  the  rate  of 
drainage  following  irrigation  on  the  determination  of  consumptive 
use.  He  assumed  that  the  drainage  rate  from  a  given  depth  of  a 
soil  is  a  function  only  of  water  content  and  independent  of  the 
rate  of  extraction  of  water  by  roots.  He  underestimated  the 
consumptive  use  when  the  drainage  rates  from  a  covered  plot  was 
deducted  from  the  total  loss  of  a  nearby  cropped  plot,  and  over 
estimated  the  consumptive  use  when  the  total  moisture  lost  from  a 
cropped  plot  was  used  as  estimate  of  the  consumptive  use.  The 
theory  of  the  instantaneous  profile  method,  referred  to  hereafter 
as  the  IPM,  was  modified  by  Rose  et  al.,  (1967)  to  permit  the 
calculation  of  water  withdrawal  by  plant  roots  as  a  function  of 
depth  and  time  with  out  neglecting  water  movement  in  the  soil. 
Using  this  method  to  separate  between  water  redistribution  and 
uptake      by      plant      roots      it      requires: 

1.  Knowledge      of      change      in      water      storage      which      is      given      by 

fZfTa 

(de/dT) dzdT 

JO    T, 

2.  An      estimate      of      the      vertical      flux      q      which      is      given      by      the 
Darcy       equation 

3.         An       estimate       of       evaporation       which       is       distinct       from 

transpiration 


T2 

E     dT     : 

Tl 


T2 
q    dT 
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where     q0     is     the     upward     flow     at     the     soil     surface     which 
occured        in        a        liquid        phase        and        that        the        vapor        flow        is 
small       in       comparison. 

4.  Amount  of  water  which  is  withdrawal  by  roots  from  a  soil 
with  an  upper  boundary  at  the  soil  surface,  a  lower 
boundary     at     depth     z     is      given     by 


r (z) dz 


0 


rT2 
(I -q-E)     dT    - 

Ti 


T2 


(ae/aT) dzdT 


T 


where  l  =  rate  of  irrigation  including  precipitation,  and 
r(z)  =  time-averaged  rate  of  water  withdrawal  by  roots  over  the 
period      T1      to      T2. 


They  concluded  that  the  method  can  be  used  succesc  fully  to 
determine  in  situ  the  pattern  of  water  «  ;  thdrawa  i  from  the  soil  by 
a      growing      crop. 

The  instantaneous  profile  technique  may  also  be  applied  to 
laboratory  flow  columns  and  has  the  advantage  of  comparing  the 
flux  in  the  still  saturated  zone  at  selected  profile  times  with 
the  outflow  rate  per  unit  area  of  the  same  times  as  measured  from 
the  volume  outflow  at  the  base  of  the  column.  Also  the 
arrangement  of  measuring  the  soil  water  content  and  soil  water 
suction  with  the  gamma  rays  absorption  and  tensiometer  pressure 
tranceducer  respectively,  permits  the  simultaneous  nondestructive 
and  rapid  response  measurement  of  soil  water  content  and  soil 
water  suction  and  provide  a  record  of  these  changes  against  time. 
Watson       (1966)       proposed       this       technique       and       applied       it       to       an 
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initially  saturated  sand  column  where  measurements  of  soil  water 
content  and  pressure  head  as  a  function  of  depth  and  time  were 
made.  In        his         investigation,         Watson        (1966)         used        the        one 

dimensional  continuity  equation  for  unsaturated  materials  to 
calculate  the  flux  at  given  time  intervals  by  integrating 
graphically  60/dT  with  respect  to  depth,  then  the  gradient  was 
determined  by  differentiating  graphically  the  total  potential 
with        respect        to        depth        at        the        same        time        intervals.  The 

instantaneous  hydraulic  conductivity  was  then  calculated  by 
dividing       the       flux       by      the      gradient. 


Before  1967,  a  problem  associated  with  the  studies  of 
separating  the  drainage  and  consumptive  use  (root  absorption  plus 
evaporation)  for  an  actual  field  situation  was  that  in  none  of 
these  studies  was  there  an  independent  measure  of  evaporation 
rate.  Van  Bavel  et  al.  (1968a)  determined  the  evaporation  losses 
from  three  precision  weighing  lysimeters  in  which  the  surface  and 
environment  were  closely  identical  to  those  of  the  test  plots. 
They  used  the  IPM  (Watson  1966)  to  calculate  the  flux  and  the 
hydraulic  conductivity  from  in  situ  measured  water  content  and 
pressure  potential  by  determining  dO/dT  and  dH/dz  at  a  selected 
group  of  time  intervals  using  a  graphical  technique  and 
integrating  de/dT  for  depth  increments  under  consideration  to 
get  the      flux,      which      was      divided      by      the      hydraulic      gradient      to 

get       the       hydraulic       conductivity       at       a       given       depth.       The       rate       of 
root       extraction,       r(z),       was       then       calculated        using 
r(z)        =        ae/dT       -        dq/dz 
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Using  this  method,  van  Bavel  et  al.  found  that  the  calculated 
root  extraction  rates  agreed  reasonably  with  the  independent 
lysimetric  measurements  of  the  water  loss  from  the  surface  to  the 
atmosphere. 

In  bare  soils,  infiltration  .evaporation,  and  deep 
percolation  depend,  in  some  measure,  upon  the  water  conent  of 
soil  profile.  Black  et  al.  (1969)  showed  that  the  cumulative 
evaporation  for  a  bare  Plainfield  sand  at  any  stages  was 
proportional  to  the  square  root  of  time  following  each  heavy 
rainfall  and  can  be  calculated  from  the  diffusivity  measurements, 
and  that  the  drainage  rate  was  an  exponential  function  of  total 
water  stored  above  a  given  depth  of  soil  profile.  Miller  et  al. 
(1971)  studied  the  effect  of  evaporation  rate  on  drainage  losses 
at  different  depth.  They  found  that  drainage  losses  increased  and 
extended  over  longer  times  as  the  evaporation  decreased,  and  that 
the  drainage  losses  at  120  cm  were  greater  and  extended  over  a 
longer     period     of     time     than     that     at     70     cm     depth     . 

Arya       et       al.  (1975)       described       another       field       method       to 

determine  K(9)  and  0(h)  which  requires,  as  well  as  the  IPM, 
either  direct  measurements  of  both  soil  water  content  and  soil 
water  pressure  head  profiles  or  the  direct  measurements  of  either 
one  of  these  variables  and  the  indirect  determination  of  the 
other  variable  through  the  separately  determined  water  retention 
curve.  In  this  method  the  boundary  between  the  upward  and 
downward  movement  of  water,  the  plane  of  zero  flux,  was 
positioned      as      it      moved      down      the      soil      profile      with      evaporation      and 
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drainage        occurring        simultaneously.  Arya        et        al.  graphically 

evaluted  the  hydraulic  head  gradient  from  total  potential, 
measured  in  the  field  using  tensiometers,  plotted  as  a  function 
of  depth  at  various  times.  Then  the  graphically  determined  values 
of  the  hydraulic  head  gradient  were  plotted  against  time  to 
determine  the  position  of  zero  flux  plane.  This  method  has  the 
same  basis  of  calculating  the  flux  as  that  suggested  by  Richards 
(1956)  and  the  same  basis  of  inferring  the  hydraulic  head 
gradient      as      that      suggested      by      Watson      (1966). 

Olsson  and  Rose  (1978)  used  the  IPM  to  determine  the 
hydraulic  conductivity  characteristics  of  a  soil  profile  that  is 
subjected  to  volume  changes  with  changes  in  water  content, 
swelling  soils,  from  in  situ  measurements  of  water  content  and 
suction  during  the  redistribution  of  water  through  the  profile  of 
a  red-brown  earth  which  exhibits  swelling  properties.  They  noted 
that  at  a  given  water  suction,  the  hydraulic  conductivity  was 
generally  lower  in  subsoil  where  micropores  dominate  these  layers 
imposing  a  high  resistance  on  water  flow  and  thus  reducing  the 
bulk      velocity      of      water      for      a      given      potential      gradient. 

Reliable  estimates  of  unsaturated  hydraulic  conductivity  are 
difficult  to  obtain  due  to  extensive  variability  in  the  field  and 
cost.  For  this  reason  equations  have  been  developed  by  Childs  and 
Collis-George  (1950),  Marshal  (1958),  and  Millington  and  Quirk 
(1959,  1960,  1961)  to  calculate  the  unsaturated  conductivity  from 
pore  size  distribution  which  can  be  characterized  easily  by  the 
standard       measurements       of       water       content       vs.         pressure.         These 
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methods  have  been  tested  with  some  success  by  JacKson  et  a). 
(1965),        Kunze        et        al.  (1968),        Green        and        Corey        (1971),        and 

Campbell        (1974).  In        developing        these        equation        it        is        assumed 

that   the  soil  is  isotropic        with        respect        to        hydraulic 

conductivity,  an  assumption  of  uncertain  validity  in  field 
soils.  van        Genuchten        (1978,        1980)        derived        a        closed        form 

analytical  solution  based  on  both  Burdine  theory  (1953)  and 
Mualem  theory  (1976a)  for  predicting  the  hydraulic  conductivity 
from  the  soil  moisture  curve.  This  model  contains  two  or  three 
independent  parameters  which  may  be  obtained  from  soil  moisture 
retention  data  by  a  non-linear  least-  squares  curve-fitting 
method  with  the  aid  of  digital  computers.  van  Genuchten  showed 
that  the  soil  water  content,  9,  as  a  function  of  pressure  head, 
h,      is      given      by: 


e    =    er    + 


©s  -   ©r 


[        1  +  (ah)n]m 


where     a,     n,     and     m     are     unknown     parameters     that     may     be     determined 
from      0(h)      curve,      and      9r      and      9S      are      the      residual      and      saturated 
value      of      9      which      can      be      determined      experimentally.        0r      may      be 
considered      as      the      volumetric      water      content      at      a      large      value      of 
pressure      head      (assuming      h      is      inversely      proportional      to      9)      where 
dH/d9      z      0      which      on      this      basis      assigns      9r     a      positive      value      and      a 
physical       basis       when       either       Burdine       or       Mualem       three       parameters 
model       is       considered.       However,       Ward       et       al.       (1983)       indicated       that 
9r     should      be      considered     as      a      fitting      parameter      rather      than      a      soil 
property,      which      allows      9r     to      be      less      than      zero.      Also      Stephens      et 
al.       (1985)       found       that       the       lower       limit       of       the       95/       confidenc 
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intervals        of        6r        was        less        than        zero,        which        is        physically 
impossible.       van       Genuchten       model       was       tested       by       van       Genuchten 
(1980),       Dane       (1980),       Ward       et       al.  (1983),       and       Stephens       et       al. 

(1985)  and  found  to  be  convenient  and  sufficiently  accurate  for 
field  application,  but  the  reliability  of  the  model  may  depend 
upon     a     reasonable     estimate     of     9r     and     95. 

The  IPM  is  perhaps  the  most  reliable  method  for  determining 
the  saturated  as  well  as  the  unsaturated  hydraulic  conductivities 
for  field  conditions.  Basis  of  this  method  is  the  Darcian-based 
flow  theory.  However,  Application  of  this  theory  to  field  is  time 
consuming,  expensive  to  characterize  the  water  storage  and 
conductivity  functions,  and  the  uncertainity  of  representation  of 
measurements  to  large  areas  owing  to  the  inherent  spatial 
variability  of  a  field  soil  (Rogowski,  1972).  Technical  problems 
and  soil  profile  charactristics  may  limit  the  application  of  the 
theory      of      the      IPM      (Baker      et      al.       1974). 

Field      hydraulic      properties      have      been      measured,      by      Black      et 
al.  (1969);        Davidson        et        al.         (1969);        Nielsen        et        al.  (1973); 

Luxmoore  et  al.  (1981);  Libardi  et  al.  (1980);  Chong  et  al. 
(1981);  Sisson  et  al.  (1980);  and  Jones  et  al.  (1984)  assuming 
that  the  moisture  content  and  suction  over  a  substantial  length 
are  uniform  which  implies  that  the  potential  gradient  is  only 
gravitational.  In  these  methods  the  soil  is  wetted  deeply  and 
allowed  to  drain  while  evaporation  is  prevented.  Under  this 
assumption,  the  method  is  known  in  the  literature  as  the  the  unit 
gradient      method.        Using      this      method      for      covered      and      non-vegetated 
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uniform       or  weakly       layered       soil       proflie       with  a       shallow       water 

table,       where       the       hydraulic       gradient       is       nearly  unity       and       water 

content     is     a  function     of     time     and     independent     of  depth,     it     can     be 

shown       that  the       average       hydraulic       conductivity,       K,       for       a       soil 

layer     between     the     soil     surface     (z0     =0)     and     depth  of     interest      (D  =  - 

z$      is      given  by 

K     =  D     *     (ee/dT) 

where  9  is  the  average  volumetric  water  content  above  the  depth  D 
and  K  is  the  hydraulic  conductivity  evaluated  at  the  water 
content  9  and  depth  D.  One  of  the  advantages  of  this  method  is 
that  the  soil  water  content  profiles,  measured  as  a  function  of 
time  following  steady  state  infiltration  conditions  in  a  field 
soil,      provide      the      only      data      necessary      to      estimate      K(9). 

Libardi       et       at.        (1980)       developed       a      method       to      obtain       K(9) 
only  assuming  a         unit  hydraulic  gradient,  the         hydraulic 

conductivity  is  an  exponential  function  of  soil  water  content, 
and  that  the  average  soil  water  content  between  soil  surface  and 
depth  L  is  a  function  of  soil  water  content.  This  method  was 
developed  by  integrating  Richards  equation  similar  to  that 
integration  done  by  Nielsen  et  al.  (1973).  They  concluded  that 
the  calculated  K(9)  may  differ  from  reality  at  each  site 
especially  with  soils  having  distinct  layers  of  horizons  of 
greatly  differing  hydraulic  conductivity.  The  limitations  of  this 
method  are  only  K(9)  information  is  obtained  and  the  necessity  to 
assume  unit  gradient  (Dane  et  al.,  1983).  The  calculated  K(9) 
values      using      Libardi      method      seemed      to      correspond      fairly      well      with 
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the  corresponding  K(9)  from  IPM  at  intermediate  values  of  9,  but 
it  did  not  correspond  at  all  at  higher  values  of  9  (Dane,  1980). 
This  behaviour  may  be  explained  by  the  fact  that  during  the  early 
drainage  period  (which  corresponds  to  low  absolute  value  of 
pressure  potential,  h,  and  higher  9  values),  at-l/dz  <0,  but  after 
several  days  of  drainage  (which  corresponds  to  higher  absolute 
value  of  h  and  lower  9  values),  aH/dz  >0  (Chong  et  al.  1981). 
However,      when      K(9)      evaluated      using      Libardi  method      was      compared 

with  K(9)  detrmined  using  the  IPM  ,  Schuh  et  al.  (1984)  noted 
that      Libardi  method      worked       well      on      coarse      and      fine      textured, 

homogeneous  materials  underlying  stratified  soil  materials,  and 
adequate  slopes  of  K(9)  curves  were  obtained  within  layered  soil 
materials  but  the  calculated  K(9)  as  a  matching  values  were  often 
inadequate.-:  suggesting  that  the  fit  between  calculated  and 
measured  values  of  hydraulic  conductivity  could  be  improved  by 
using        field-saturated        K(9)        values. 

Jones        et        al.         (1984)        compared       five        different        methods, 
including        Libardi  method,        to        estimate        K(9)        from        in        situ 

measurements  assuming  a  unit  hydraulic  gradient.  They  concluded 
from  their  study  that  the  five  approximate  methods  were  useful  in 
developing  a  rapid  and  rough  estimate  of  soil  water  properties 
over  large  areas,  but  the  methods  were  not  as  useful  for  a 
particular  location  where  soil  water  properties  need  to  be 
precisely        known. 

Ahuja      et      al.       (1980)      proposed      a      method      to      determine      the 
hydraulic       conductivity       as       well       as       the       soil       water       characteristic 
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by  analyzing  of  the  drainage  phase  tensometric  data,  combined 
with  the  field  measured  value  of  near  saturated  hydraulic 
conductivity  and  one  soil  moisture  sampling  during  drainage.  This 
method,  when  applied  for  data  obtained  at  five  sites  chosen  for 
maximum  variation  of  materials  between  and  within  sites,  resulted 
in  a  good  estimate  of  K(h)  as  determined  with  the  IPM  (Schuh  et 
al.         1984).  Unless        an        accurate        determination        of        saturated 

hydraulic  conductivity,  from  steady  state  infiltration  prior  to 
drainage,  is  used  as  a  matching  value  as  suggested  by  Ahuja  et 
al.  (1980)  the  matching  value  should  be  representative  of  the 
actual  log  K  vs.  log  h  (where  h  is  the  pressure  potential)  to 
avoid       inaccurate       results       (Schuh       et       al.        1984). 
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ABSTRACT 

Three  K(6)  functions  were  used  to  describe  hydraulic  conductivity 
data  from  a  layered  field  soil.  The  functions  were  <o(e/em^- 
Ko((e"ec)/(©m-©c))n  and  K0exp{a(e-em)i  where  K0  =  K0(z)  =  g(z)Km, 
g(z)  is  a  scaling  factor  that  varies  with  depth  z,  Km  a  constant 
and  ec,  em,  0,  n,  and  a  are  parameters.  For  each  function  three 
cases  were  considered:  Case  1  fit  discrete  values  by  depth  to 
g(z)  and  a,  &  and  n:  The  second  case  treated  the  scaling  and 
exponential  parameters  as  continous  functions  of  depth:  And  the 
third  case  fit  scaling  factor  as  a  discrete  function  of  depth  but 
held  0,  n  and  a  constant.  More  variation  in  r2  and  MSE  was  found 
between      cases      than      between      functions. 


